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The internal structure of hadrons is important for a variety of topics, including the hadron form 
factors, proton spin and spin asymmetry in polarized proton scattering. For a systematic study 
generalized parton distributions (GPDs) encode important information on hadron structure in the 
entire impact parameter space. We report on a computation of nucleon GPDs based on simula- 
tions with two dynamical non-perturbatively improved Wilson quarks with pion masses down to 
350MeV. We present results for the total angular momentum of quarks with chiral extrapolation 
based on covariant baryon chiral perturbation theory. 
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1. Introduction 

The internal structure of nucleons has attr acted much attention in the contexts of the nucleon 
form factors, proton spin and spin/charge asymmetry in deeply virtual Compton scattering and 
so on. For a systematic study of the nucleon internal structure, generalized parton distributions 
(GPDs) are introduced through the off-forward matrix elements of quark-bilinear operators: 

j^e^^^{PU-^)y'^qmP)=mP'){y'H{x,^,t) + i^E{x,^ (1.1) 

with a light cone vector n and the momentum transfer A = — P as functions of the quark momen- 
tum fraction x, the skewedness t, = — « • A/2 and the virtuality t = A^. The axial counterparts are 
denoted by H and E. Since the GPD is defined with the finite momentum transfer in contrast to the 
conventional parton distribution functions, partons bring us the informations on hadron structure in 
the transverse space. 

In this contribution, we report on the first moments of GPD, so called generalized form factors, 
for nucleon, as a function of the virtuality calculated on the lattice with unquenched configurations 
of QCDSF/UKQCD collaboration. 

In the forward limit these generalized form factors provide the total angular momentum of 
quark in the nucleon through Ji's sum rule [|l|], 

f = ^ j'^ dxx{H (x, § , 0) + £ ^ , 0) ) = ^ (A20 = 0) + B20 (? = 0) ) . (1.2) 

Combined with the quark spin contributions to the nucleon obtained as the forward value of the 
axial form factor, 

= \j dxH{x,^,0) = h,o{t = 0), (1.3) 

we compute the orbital angular momentum of quarks as = 7^ — 5^. Using the results of chiral 
perturbation theory (;tPT) for chiral extrapolation to the physical point, we discuss the angular 
momentum carried by quark in the nucleon. 



2. Generalized form factors on the lattice 



The Mellin moments of the GPDs are known to be expressed by polynomials in terms of ^ 



dxx' 



n-l 



Hix,^,t) 
E{x,^,t) 



[in-l)/2] 
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k=Q 



Ik 



An,2k{t) 
Bn,2k{t) 
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(2.1) 



The generalized form factors An^2k,Bn^2k and C„ are defined from the coefficients of this expansion. 
Since the integration by x makes the quark operator local, the {n — l)-th moments can be calculated 
on the lattice through the matrix element of {P'lqy^^^D^^ ■ ■ ■D^"^q\P) by taking a ratio of the 
three- and two-point functions. 

To estimate these correlation functions, 400 to 2200 configurations are used for each j8 , K with 
two flavor Wilson fermion with the clover improvement. Simulations are performed with various 
set of parameters j8 and K corresponding to the lattice spacing less than 0.09fm and pion mass cov- 
ering from order of IGeV down to 350MeV with a reference scale ro = 0.467fm. Nonperturbative 
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Figure 1: The isoscalar axial form factor Aio(f) for j3 — 5.29, K — 0.13632 with dipole fit (left) and the 
forward values with jPT fit (right). The open star in the right panel represents the latest experimental value 
of HERMES. 



renormalizations are incorporated to convert the lattice results into the values in the MS scheme at 
a scale of pL^ = 4GeV^. The ff{a) improvement of the quark energy-momentum tensor are carried 
out through the boosted perturbation theory following ref.[Q] and the tadpole improved version 
is used for the axial current following ref.[^. We note that the contributions from disconnected 
diagrams are not included in the present lattice results. 



3. Lattice simulation results and chiral extrapolation 

We focus on the generalized form factors A20 and B20 as well as the axial form factor A\q to 
evaluate the quark angular momentum in the nucleon. 

Typical t dependence of the axial form factor in the isoscalar channel are shown in Fig|I[ The 
obtained lattice data agrees well with a fitting by the dipole form, Aio(O) /(I —t/rn^Y. The forward 
values obtained by setting ? = present a smooth pion-mass dependence as shown in the right panel 
ofFig.[T[ 

Here we use an expression derived in a heavy baryon ;^PT |^, 



Hd,n^ , 3m|gl ( ml \ 



A;:+^(0) = a„,, 



-Pn,kml + ir{ml), (3.1) 



for the chiral extrapolation with fitting parameters aio and j3io at a scale of A = IGeV. As the heavy 
baryon formalism is valid only for the small pion mass, we restrict the data points at pion masses 
less than 500MeV. Then it turns out that the chiral log term gives a strong mj^ dependence for small 
rrijc region and the extrapolated value is eventually comparable with the latest experimental value 
of deep-inelastic scattering reported by HERMES [^]. With this extrapolation, we obtain the quark 
spin contribution in the nucleon as A^q '^(0) = 25'"+'' = 0.402 it 0.024 at the physical pion mass. 

Next we show a typical t dependence of generalized form factors A2o,B20 and C2 in the 
isovector channel in Fig.^. Up to IGeV^, the lattice results of A20 agree with the dipole fit 
A20 (0) / ( 1 - ? /ml f. The dipole mass niD of A20 shows a smooth mji dependence shown in the 
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Figure 2: Generalized form factors in the isovector channel for j3 = 5.29, K = 0.13632 with dipole fit for 
A20 and the dipole mass of A2q- The open star represents the experimental value of /2 tensor meson mass. 



right panel of Fig.Q We see that the dipole mass seems to extrapolate to the observed mass of 
tensor meson /2 at the physical point. This fact contrasts with a mass scale of the electromagnetic 
form factors comparable with the vector meson mass [||]. 

However, the empirical dipole fit for the generalized form factors has no solid justification 
from a theoretical point of view. Therefore we count on the covariant baryon chiral perturbation 
to extract the forward values of B20 as well as the chiral extrapolation of A20 and 620- 

The forward values of A20 are identical to the quark momentum fraction (x)^ and are shown 
as a function of m| in Figj3[ Both in the isoscalar and isovector channel, the lattice results show 
a moderate pion mass dependence. These values are extrapolated to the physical point by the 
following eqs. derived in the baryon xVT, 
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Figure 3: 

open stars 



The forward values of A20 in the isoscalar (left) and isovector (right) channel with fits. The 
represent the phenomenological values from CTEQ6. 
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Figure 4: The forward values of B20 extrapolated by jPT in the isoscalar (left) and isovector (right) channel 
with xV\ fits. 
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for the isovector channel, where Mq is the nucleon mass in the chiral limit. We perform 2-parameter 
fits with aIq, C9 for the isoscalar and a^^, cs for the isovector channel at a scale of A = IGeV and 
fixed the other values following ref.[^. A strong m^t dependence is observed especially for the 
isovector channel, but the extrapolated values in both channels overshoot beyond the phenomeno- 
logical values estimated using the CTEQ6 parton distribution functions. The chiral extrapolation 
gives A^+'* = {xY+'^ = 0.572 ±0.012 for the isoscalar channel andA^^'* = {x)"-'^ = 0.198 ±0.008 
for the isovector channel at the physical pion mass. See ref.|10] for discretization effects of these 
form factors. 

In contrast to A20, the forward values of S20 cannot be calculated directly from the lattice 
simulation since the kinematic pre-factor for B20 vanishes at zero momentum transfer. Again we 
use the expressions of covariant baryon ;t;PT [^], 
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to fit the lattice results as a function of t and m^r. Free fit parameters are Z?2o > ^rid 8^)' for the 
isoscalar and isovector channel respectively, and the parameters of d^Q fixed by the fitting of 
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Figure 5: Total angular momentum of quark in the nucleon with jPT fit (left) and spin, orbital angular 
momentum of quarks (right). The open symbols represent the extrapolated values to the physical pion mass. 



A"q^. The forward values of B20 extracted from this fit with fixed m,^ are shown in Fig.^, where the 
solid lines represent the section of fitting surfaces at ? = 0. 

Since the forward value of is equivalent to the difference of 27"? — (x)^, the small values 
of the lattice results indicate the cancellation between total angular momentum and momentum 
fraction of quarks. However, the ;^PT fit suggests a sizeable bending through the chiral extrapo- 
lation, which makes a sharp contrast to the chiral quark soliton model [pi]]. We obtain '^(0) = 
-0. 120 ± 0.023 at the physical point for the isoscalar channel and B20 '^(0) = 0-269 ± 0.020 for the 
isovector channel. 

Combining all these data, we can estimate the angular momentum of quarks in the nucleon. 
The pion mass dependence of the total angular momentum 7^ for u and d quark is shown in Fig.§. 
The strong dependences of A20 in the isovector channel and B20 in the isoscalar channel are en- 
hanced for the u quark and gives a significant suppression of J" near the physical point, while these 
dependences cancel each other for the d quark and so keep the value of J'^ small. The extrapolation 
to the physical point gives r = 0.230 ± 0.008 and J'^ = -0.004 ± 0.008. 

The quark spin s'^, the total and orbital angular momentum U = J'^ — are shown in the right 
panel of Fig.|5[ We obtain the values at physical point by the chiral extrapolation as 7"+'* = 0.226 ± 
0.013, = 0.201 ±0.024 and L"+'* = 0.025 ±0.027. The results show that the total angular 
momentum of quark is comparable with the quark spin and hence the orbital angular momentum is 



consistent with zero, which agrees with the result of ref.[ 12]. 



4. Conclusions 

We have carried out lattice simulations to calculate the first moments of GPDs, which play an 
important role for the proton spin, quark transverse density and deeply virtual Compton scattering. 
The lattice results of the generalized form factor A20 are fitted to the dipole form for small —t region 
and the dipole mass turns out to be comparable with the observed tensor meson mass. Since this 
empirical fit has no solid justification from a theoretical point of view, we employ baryon ;^PT to 
take the forward limit of B20 and to chirally extrapolate the form factors to the physical pion mass. 
The resulting values indicate that the total angular momentum of quarks in the nucleon is of the 
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same size as the quark spin contribution, while the orbital angular momentum is consistent with 
zero. 



Further analyses are needed to estimate the finite size effects ||13|], contributions from discon- 
nected diagrams and so on. Results with lighter pion masses will be forthcoming. 
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